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Once the quark chemical potential fi is introduced in finite density QCD, the fermion determinant appeared 
in the path integral measure becomes complex. In order to investigate the phase effect of SU(3) lattice QCD 
(2-flavors), we calculated the fluctuation of the phase of det A(/x) on a 8 3 x 4 lattice at /i = 0.1 and 0.2. Then we 
calculated the chiral condensate and Polyakov line in the no phase and reweighted cases. There is little difference 
between these two cases at fi — 0.1 and 0.2. We consider a possible reason for this result below fi < 0.27 in terms 
of Zz symmetry. 



1. Introduction 

Finite density QCD £Q has attracted consid- 
erable attention in high energy physics, nuclear 
physics and astrophysics. Many theorists now be- 
lieve that QCD has a very rich structure when we 
study it in temperature and density parameter 
space, and some experimentalists want to reveal 
it. 

Fodor and Kaiz 2 studied QCD at T ^ 
and at ^ by the multiparameter reweighting 
method. If two measures at fi = and /j, ^ over- 
lap well, then the method works well. However, 
for large /i the method may not work due to the 
phase effect. Better overlap may be obtained by 
a simulation with a phase-quenching measure[3]- 
Here the difference between two measures is a 
complex phase e %e . If the phase fluctuates very 
much, then the method again fails, which is ex- 
pected to happen for large [i. 

It is therefore important to investigate the be- 
havior of the complex phase in (T, fi) parameter 
space. In this study, we calculate the complex 
phase for various and fj, and study the phase 
fluctuation. Furthermore, we use the complex 
phase to obtain results at finite 
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2. Numerical Investigation 

The lattice partition function for KS fermions 
is given by 

Z = jDU{detA(n)) Nf/i e- fJS °. (1) 

We use the phase-quenching measure ~ DU 
| det A(/i)| JV - f / 4 e~' 5Ss in our Monte Carlo simula- 
tions. We perform calculations on a 8 3 x 4 lattice 
at to = 0.05 for Nf = 2 using the R-algorithm. 
The trajectory length is set at 0.5 with a step 
size of At = 0.01. The first 1000 trajectories 
were discarded for thermalization. 

For Nf = 2, the expectation value of an oper- 
ator O is given by 

(O) = ^ J DU(det A(/i)) 1 / 2 0e-' 3S » 

J DU\ det A^e^^Oe^ 
~ JDU\ det A| 1 /2 e ie/2 e -0S g 

_ JPUldetA^e^Oe-^" 
JwjdetAfVV^ 

JPU\det A|l/2 e »9/2 e -/35 9 _ ( 0e »9/2) Q 

/ J DU\ det A^^e-^^ ~ (e^/ 2 ) ' (2) 

where (. . .)o is the expectation value with the 
phase-quenching measure. 
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Figure 1. (cos(£?/2)} f° r A* 
functions of j3 



0.1 and 0.2, as 



Figure 2. No phase and reweighted chiral conden- 
sates (Re-0V>) for [a = 0.1 and 0.2, as functions of 



To obtain the phase 9 of the fermion deter- 
minant, it is necessary to calculate the fermion 
determinant. We consider two methods for cal- 
culating det A (fi) : (I) calculation of eigenvalues, 
(II) LU factorization. We found that LU factor- 
ization is much faster than method (I). Thus we 
have decided to use method (II) for the phase 
calculations. Figure shows the fluctuations of 
(cos(0/2))rj. The fluctuations increase with in- 
creasing [i. The fluctuations decrease for [3 > (3 C , 
where critical coupling (3 C is around 5.3. 

Using the phase 8 of detA(/i), we estimate 
the reweighted values of the chiral condensate 
(Rcipip) according to reweighting formula J3J), 



i.e., 



(Re inp) 



(Reipip cos(0/2)) o 



(3) 



(cos(0/2)) o ' 

Figure |21 shows the graphs of the no phase 
case and the reweighted case. Values for the 
reweighted case are very similar to those for the 
no phase case. The latter agrees with the /i be- 



havior reported by Kogut and Sinclair 0|. 

Figure [3] shows Polyakov lines for the no phase 
case and the reweighted case. The two are also 
almost the same. In our approach, the conjugate 
gradient method does not converge at fi > 0.27 
and [3 = 5.20. We are performing simulations 
below fi c . If we can perform simulations above 
/i c , we might see the difference between the two 
cases. 



3. Z% symmetry 

Why does the phase effect not appear? Aarts et 
al. |S] proposed a partial summation in configura- 
tion {{/} connected by Z3 symmetry. Let us write 
a set of link variables {ZU} after the Z% transfor- 
mation, i.e., we multiply z — exp(2iit/i) S Z% = 
{l,e 2i7r / 3 ,e 4 "/ 3 } to all {U 4 (x)} on a time slice. 
Since z is an element of SU(3), D(ZU) = DU. 
The gauge action is invariant, Sa(ZU) = Sq(U). 
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It is interesting to see whether the phase disap- 
pears in E. We are now trying to calculate this 
quantity. 

4. Summary 

We calculated the fluctuation of the phase of 
det A(/x) on a 8 3 x 4 lattice at fi — 0.1 and 0.2. In 
the region of small f3, the fluctuations are large. 
We calculated the chiral condensate (Re ijji/j) and 
the Polyakov line (ReL) in both no phase and 
reweighted cases. There is no difference between 
them. Wc are currently running at fi = 0.25, 
varying (3 through the crossover region. 

We are trying to recognize these no phase ef- 
fects below // < 0.27 in terms of Z% symme- 
try. For more larger it is expected that this 
reweighting method will be effective. We must 
try to calculate the case of /i > 0.27. The calcula- 
tions of gluon energy density and fermion number 
density are also in progress now. 



Figure 3. No phase and reweighted Polyakov lines 
(Re L) for /i = 0.1 and 0.2, as functions of /3 



Therefore, 

Z = J D(ZU){det A(ZU)) 1/2 e- ps <> {zu) 

W(det A(Z[/)) 1 / 2 e~ /3 ' s ^. (4) 



The same is true for {Z 2 U}. Then we obtain 
Z= X - I DU (det A(U)) 1/2 



( /det A(ZU)\i/2 / 

( 1 + (imWT) + ( 



det A(Z 2 C/)^V2 
det A (IT) 



-ps g (U) 



(5) 

For any configuration {U}, {ZU} and {Z 2 U} will 
appear in a whole set of ensembles. Therefore, it 
is worth calculating the following partial summa- 
tion in configuration {U} connected by Z3 sym- 
metry: 

/det A(Z£T)\i/2 /det A(Z 2 C/)\V 2 
= + V detA(£7) ) + \ detA(C/) ) 
= \E\e ie . (6) 
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